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We investigate the newly discovered supersolid phase by solving in random phase approximation 
the anisotropic Heisenberg model of the hard-core boson 4 He lattice. We include nearest and next- 
nearest neighbor interactions and calculate exactly all pair correlation functions in a cumulant 
expansion scheme. We describe the properties of the normal solid and supersolid phases and argue 
that based on analogies to the fermionic half filled extended Hubbard model the supersolid state 
corresponds to a bond-ordered-wave. 
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Owing to its theoretical and experimental importance, 
the observation |l| of non-classical rotational inertia asso- 
ciated with co-existing solid and superfluid (SF) phases 
generated a lot of experimental controversy 0, IH and 
theoretical debates 0, H, However, the latest mea- 
surements 0] clearly show that there is an actual super- 
solid (SS) phase. Ref. poses serious limitations for 
previous theories [4|, [5( but simultaneously open the door 
for systematic theoretical approaches, which is the goal 
of the present Letter. 

Since the work of Penrose and Onsager Q it has been 
known that the SS has to exhibit simultaneously two 
types of order, namely diagonal long-range order (DLRO) 
associated with the periodic modulation in a crystal and 
the off-diagonal long-range order (ODLRO) associated 
with the phase order in the condensate. Well-known 
for fermionic systems such as excitonic insulators [8|, 
but, for bosons the difficult reconciliation of such puz- 
zling behavior has been debated since the late 1960's 
0,[lO|. After Kim and Chan's [l[ discovery however, only 
two kinds of theoretical approaches have been pursued, 
namely numerical simulations 0] and phenomenological 
descriptions [B|, [||, while many-body theories, following 
the pioneering works of Matsuda and Tsuneto; Liu and 
Fisher [icj |. have been completely neglected. 

In this Letter we fill this gap, presenting a theory be- 
yond the mean-field of SS. We use a cumulant expansion 
and solve the RPA equations for the Green's functions, 
taking into account the pair correlations rigorously. We 
will clarify the controversy over the role of the vacan- 
cies and defects, which have long been proposed to have 
a crucial role in the formation of a SS phase. We find 
that vacancies and interstitials are present even at zero 
temperature in the SS phase, both condense and SS may 
be regarded as a bond-ordered-wave as it exhibits alter- 
nating strength of the expectation value of the kinetic 
energy term on bonds. Also, our model confirms that 
the SF to SS transition is triggered by a collapsing ro- 
ton minimum; however, the SS phase is stable against 
spontaneously induced superflow. 

The quantum lattice gas (QLG) model as given by 



Matsubara and Matsuda [l2j is K = Hqlq — fiN, with 

K = fl^2 Tli + ^ U ij( a l ~ a ])i a i ~ Qj) + VijUiUj , 

(1) 

where the indices i and j run over all lattice sites, j are 
the nearest and next nearest neighbor hopping param- 
eters and Vij takes nearest and next nearest two parti- 
cle interactions into account. The creation and annihila- 
tion operators a] and represent hard core bosons, i.e., 
4 He-atoms. Originally, Eq. Q] was used for SF [l| and 
later to study the possible coexistence of DLRO and 
ODLRO for SS. We start from a bcc lattice, which we 
separate into two sub-lattices, see, Fig. [T]a), establish- 
ing a natural way to define DLRO of solids: sub-lattice 
A represents the 4 He ions, while sub-lattice B the in- 
terstitial centers. A normal solid (NS) then is given by 
a fully occupied sub-lattice A, while in a liquid phase 
the occupation number on both A and B sub-lattices are 
equal. 

Conventional techniques of quantum field theory are 
not applicable to the QLG model due to the lack of 
Wick's theorem for hard-core bosons. Therefore we 
take advantage of the well known 1-to-l correspondence 
with spin-1/2 models imposed by the following identities: 
a] = Sf - iS], a , = .Sj + iS] and //, =l/2-S?. This 
mapping is exact and transforms the QLG model to an 
anisotropic Heisenberg model (AHM) in an external field: 

H = h z ^S* +^ J ij S t S j + X! X! ■ l 'j S '' S - ■ ( 2 ) 

where 4 = = -2 Uij , h* = 4~E, 4" 

These J's have to be specific for He. The interac- 
tions between the 4 He atoms are controlled by van der 
Waals forces and their repulsive nature at very short 
distances determines negative nearest neighbor interac- 
tion j| , evoking AF ordering in the spin language. The 
corresponding Lennard-Jones potential is short ranged 
and therefore it is sufficient [10|] to only consider near- 
est and next nearest neighbor interactions. Hence, j| = 
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FIG. 1: a) A bcc lattice consists of two interpenetrating sc 
sub-lattices, A (continuos line) and B (long-dashed line). For 
simplicity we only drawn the two dimensional case. In order 
to connect our result to the known bond-ordered-wave struc- 
tures, we reduced the bcc lattice into a sc lattice (dotted line) 
by unfolding the Brillouin zone, b) The obtained bond-order- 
wave in the [100] direction is shown. The continuous (dashed) 
lines the calculated expectation values are higher (lower). 

-qijfcAjeB? A = -Q^JieAjEA' J i = ~<li- J i£AjeB and 
jj = —q2-Jj e AjeA where q\ = 6 and q2 — 8 are the 
number of nearest and next nearest nei ghb ors on the bi- 
partite bcc lattice. For 4 He we will use [10j, = 1AK, 
Jj = 0.5K, jf = -3.8K and jj = -1.7K. The results 
do not change if the j" 's are within ±2 range of these 
values and Jg > J'i, J T 's values remain positive. Differ- 
ent values of the J's will lead more/less robust SS phase, 
i.e., the SS phase will occupy a larger/smaller area in the 
global phase diagram of 4 He. 

The ferromagnetic (FE) and anti-ferromagnetic (AF) 
phases do not exhibit any off-diagonal long range order: 
they correspond to the normal fluid and NS states of 
QLG, respectively. The remaining two phases of AHM 
are the canted ferromagnetic (CFE) and canted anti- 
ferromagnetic (CAF) states. As the spins in the canted 
phases are tilted there exists a non- vanishing transversal 
component, i.e., {cii} is non-zero. This is the usual or- 
der parameter for SF; hence CFE corresponds to SF and 
CAF corresponds to the SS phase, respectively. 

The Heisenberg model has been studied thoroughly al- 
though there exist only classical mean-field solutions for 
the CFE and CAF [loj]. We solve the equation of mo- 
tion in RPA to calculate the time-temperature-dependent 
Green's function G^(t) Ret/Adv = T»0(±t)([Sf (t), S%])], 
for all four phases. In order to preserve pair correlations 
as accurately as possible, we have chosen a cumulant 
decoupling scheme to approximate higher order Green's 
functions. This is the first time that a thorough cumu- 
lant RPA many-body calculation is solved for the QLG 
and/or AHM model and as such it represents a vital step 
in the understanding of these models and their phenom- 
ena. 

The cumulant decoupling is based on the assumption 
that the last term of the equality {ABC) = (A)(BC) 



+ (B)(AC\ + (C)(AB) -2 (A)(B)(C)+ ((A - (A))(B - 
{B))(C — {C})) is small and thus negligible [13]. By using 
this we calculate exactly all the pair correlations. This 
decoupling scheme couples six Green's functions, one for 
each spin component in x, y and z direction on the two 
sub-lattices respectively, to a set of six equations. Due 
to the enormous number of terms (1024 in total) within 
these Green's functions we do not intend to reproduce 
their exact form. It can be shown or alternatively argued 
that the Goldstone theorem of gapless modes imposes an 
additional condition [3] on the mean fields in the SS (and 
SF) phase, reducing the number of order-parameters by 
two. As this condition does not apply to the NS phase 
their Green's functions are structural different. 

In the next step we calculate the phase diagram of 
AHM and derive all thermodynamic quantities, e.g., the 
pressure, the entropy and the specific heat; which are 
of particular interest. Note that, by fixing the J's, the 
chemical potential, [i, remains the only variable in the 
model, which we use as a fitting parameter. This is a 
natural choice as fi is related to the pressure; they are 
inversely proportional. In order to derive a rigorous re- 
lationship between these two parameters we note that 
Qqlg — = Qahm, where G is an arbitrary potential, 
e.g., the Gibb's free energy. Then the Maxwell equation 
gives the desired relationship (dp/dfj,)T,v = {dN/ dV)r,n 
= Nq(1 — e)/V, where N and No denotes the number of 
particles and lattice sites, respectively. Here we take into 
account explicitly the net vacancy density, e, as it has 
been long proposed that vacancies and defects may play 
a crucial role in the formation of the SS phase. 

The net density of vacancies i.e. the number of vacan- 
cies minus the number of interstitials: e = 1 — ua — ns = 
(S A ) + (S g ), when taken as the measure of incommensu- 
ration has been recently proposed [|| as the key parame- 
ter to characterize quantum crystals. Existing numerical 
calculations could only follow individual atoms [3] and 
as such found very difficult to pick out the effects of the 
density of vacancies. But this is not the case in our many- 
body RPA approach. In this Letter we investigate for the 
first time the effects of s on quantum solids in general and 
SSs in particular. 

In order to understand the physical origin of the dif- 
ferent phases of the QLG model, we study first the quan- 
tum fluctuation at T = 0. In previous approaches [13] 
such an analysis was impossible to perform as the total 
spin magnitude always is locked in, or restricted to the 
vicinity of total spin 1/2. However, in our case the short 
and long-range correlations are taken into account by the 
cumulant expansion; we can explicitly see the quantum 
fluctuations at T = 0. In Fig. [2] we plotted the magni- 
tude of the total spin and the angle of the spin relative to 
the x - y - plane as a function of the external magnetic 
field h z . At large h z (/i, or small pressure) the AHM 
model is locked into a ferromagnetic phase; 4 He is a nor- 
mal liquid. There is no competition between S A and Sg 
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FIG. 2: The magnitude of the total spin on site A and B and 
the relative angle «a/b = tan -1 {S A / B )/{S A i B ) is plotted at 
T = in a) and b), respectively. In a) \Sa\ and \Sb\ are 
identical except for the SS (canted anti-ferromagnteic) phase. 
In b) the individual a's are plotted (long dashed line) and 
their differences (continuous line). 



and the magnitude of the total spin is always 1/2, i.e., a 
fully polarized state appears for any h z . 

As h z (n) is lowered (increasing the pressure) , the con- 
tribution of the transverse term SfSJ + SfSj, where 
i,j — A,B, increases until (S A / B ) an< ^ (&a/b/ become 
non-zero, giving rise to ODLRO. Consequently (S A , B ) ^ 
1/2. This is the well-know A transition, where a fully po- 
larized FE phase becomes canted. From Fig. [2] we can 
see that the magnitude of the spin is a smooth function of 
h z while the relative angle between the A and B spins is 
still vanishing, i.e., there is still no competition between 
the A and B sites. In the language of the QLG model, 
the emergence of the CFE phase means n A = rig ^ 
and, being at T = 0, condense. This is the well-known 
SF phase and hence we do not go into further detail but 
merely emphasize that within the AHM decsription SF 
appears as a competition between h z S A ^ B and the tran- 

verse coupling, S X A/B S A/B + S V A/B S V A/B . 

As we further lower h z (/i), i.e., as we increase the pres- 



sure, h z S A / B starts to compete with the Ising coupling, 
S z A j B S z A ^ B within the already well developed ODLRO. 

This competition leads to the CAF (SS) phase, where j" 
dominates over the effects of the external magnetic field. 
From Fig. [2] we can see that the magnitude and relative 
angle of the A and B spins differ. The relative angle 
variation gives the CAF phase, while the difference in 
their magnitude is responsible for the non-zero net den- 
sity of vacancies. n A / B = 1/2 — (S A , B ), Fig. [^implies 
that < n A < 1/2 and 1/2 < n B < 1 in the SS phase. 
Hence, as we increase the pressure, there will be a density 
transfer from the vanancies to the interstitial sites, i.e., 
n a — ► n>B- This particle density transfer is induced by 
the increasing strength of pair correlation (S A S B ) and 
it only happens in the SS phase; hence it stands as the 
key in understanding the origin of this phase. Within 
the QLG model, this effect will result in a alternating 



strength of the expectation value of the kinetic energy 
term on the bonds. This dimer order is called a bond- 
order-wave (BOW) in the half filled fermionic extended 
Hubbard models [15j. Calculating the expectation value 
of the kinetic energy of the QLG model we plot the [100] 
BOW corresponding to SS, see Fig. [T]b). Similar BOWs 
appear also in [010] and [001] corresponding to the roton 
minima (see below). Hence, we conclude that the boson 
condensation phenomenon of SS is a BOW. 

Since the 1970's it has been suggested that a SS 



transition might be triggered by a collapsing roton min- 
imum. In our approach we could verify this for the first 
time and found that the excitation spectrum indeed goes 
soft at [100] ([010] and [001]) at the SF to SS transi- 
tion. This is true for input parameters which satisfy 



J i + J 2 + J i — J 2 ' < 0. In all these cases the spectrum 
is not effected in the [111] direction. There are situa- 
tions when the roton minimum collapses at [111] instead 
of [100], nameley for -2J 2 T /(J 1 T + J 2 T + jj) > 0. How- 
ever, these cases induce structural instability and refer 
to a transition to a different state. 



It was also conjectured [Tl| recently that the SF phase 
in the vicinity of the SF to SS transition is unstable 
against spontaneously induced superflow and superflow 
associated with vortices. This can be easily shown not 
to be true in our RPA approach. A net superflow is 
either given by a moving condensate which results in 
a gradient of the phase of the wave-function or equiv- 
alently by a moving environment while the condensate 
stays at rest. Using the latter, we add to Eq. (Q]) 
the term, J d?xtp' (x)(iKv n ■ V)i/)(x) which transformed 
into the AHM language gives an additional contribution 

Ey J ij( S f S j ~ S j S I) t0 E q- ©■ Here the nearest and 
next-nearest neighbor cross coupling constants are anti- 
symmetric (Jy = — J*) and are zero for directions per- 
pendicular to the motion of the environment v n . Cal- 
culating the excitation spectrum it is noticed that even 
thought J x gives a small contribution to the spectrum, 
it does not influence the roton dip that triggers the SF 
to SS transition. 

At finite temperature we analyze the NS to SS transi- 
tion first, as this is known experimentally [l[. At constant 
chemical potential, the NS to SS transition is of second 
order with a discontinuity in the net vacancy density. 
In Fig [3] a) we show a typical curve of the net vacancy 
density at constant chemical potential. We have chosen a 
value for the chemical potenatial corresponding to a pres- 
sure just above the critical pressure measured by Kim 
and Chan |lj. The discontinuity in e is suggestive for a 
commensurate-incommensurate transition, when the net 
vacancy density is taken as a measure of incomensura- 
tion, as it is done in Ref. [6|. 

The thin dashed curves in Fig[3]a) show what the net 
vacancy density would be if a SS phase would not appear: 
at T — > the ground state is a perfect, commensurate 



4 




T/T 



FIG. 3: a) Net vacancy density at the NS to SS transition for 
fixed h z . The dashed curves are what the net vacancy density 
of the normal crystal would be without SS. Inset shows the 
T — limit as a function of h z . b) Net vacancy density 
of 4 He (left axis) for different values of pressure (right axis) 
superimposed onto Kim and Chan's [lj phase diagram (filled 
circles) . 

crystal, i.e., a Mott insulator, which is not a SF. In NS the 
density of vacancies increases faster than that of intersti- 
tials with increasing temperature, as reported in Ref . [l6| . 
indicating that thermal fluctuations favor vacancies more 
than interstitials. This is not the case for the SS phase, 
where all three e's have the same mean-field SF form: 
e (l - \J A/ksT) exp(-A/fc s T). Here, e = e(T = 0), 
shown in the inset of FigQ]a), shows that the SS phase 
is driven by the chemical potential, similarly to a corre- 
lation driven Mott transition: eq oc (/i — /i c ), where the 
proportionality constant for our model is ((S%)/{S%)). 

Recalculating the NS to SS transition as a function of 
the pressure and superimposing it onto the known phase 
diagram [lj], we obtained Fig[T]b). It shows a very good 
fit to the experimental points. For these isobar curves, 
as it can be seen, the net vacancy density is constant in 
the SS phase. This is because the thermal expansion of 
the SS phase is evanescently small. 

In conclusion, we have shown that the transition into 
a SS phase is of the incommensurate-commensurate type 
with properties suggesting a BOW. The AHM model in 
RPA does explain most of available experimental data of 
the SS 4 He phase, which we will present in a future pub- 
lication. Due to the interest in the recent experiments of 
Aoki, et al., Q, we briefly outline the solution to Ref. 2] 
based on Eq. O As we have shown, the SS of 4 He is 
a true solid hence, the underlying lattice is the reference 
frame. In torsional oscillator measurements, however, the 
reference frame is accelerated. Hence, the kinetic energy 
will acquire a new cross-coupling term, as presented ear- 
lier: J x a J2i j (SfSj — SjSf), where a is the acceleration. 
This will contribute with J T ^ . cos(0j — <f>j + J x / J T a) 
to the vortex dynamics. This implies that the mo- 
tion of a single vortex/antivortex is given by <ir v /dt oc 



±fi x dv n /dt, rather than the well-know SF equation 
e?r v /dt oc ±f2 x v n , where v„ = r x f! is the velocity 
of the rotating cylinder = f2rjCos(u;t). Thus, for a SS 
the signal p s /p is independent of the frequency as seen 
in the Aoki, et al., experiments i.e., the acceleration 
and not the rim velocity is the defining parameter in the 
motion of vortices. 
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